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SUMMARY 


A comparison is made between experiment and theory for the calcu¬ 
lation of shielding by a structure against fallout radiation. 

Experimental results for exposure penetrating a roof slab, and 
for reduction factor, are found to agree with moments theory calcu¬ 
lations, often to better than 10% when the geometry factor L (u),X) 

a 

was used. 

A comparison between experiment and theory may be inaccurate due 
to anisotropy of the experimental source, to lack of source reflection 
in roof experiments, and due to error in estimating the thickness of 
the roof. The magnitudes of the errors due to these effects are inves¬ 
tigated, and found to be small, but not necessarily negligible. 

Detailed results of penetration in iron and concrete, due to plane 
isotropic sources of various energies are given, and a review of ground 
contribution to fallout radiation penetration through simple structures 
is appended. 
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1 INTRODUCTION 


Standard methods of calculation of the protection afforded by a 
structure against fallout radiation depend upon the separation of the 
solution into a barrier and a geometry factor. These factors have been 

calculated in large part by moments theory, and have been presented as 

. ( 1 , 2 ) 

graphs 


The primary aim of the research project reported herein was to 
determine how valid the method of calculation is when applied to the case 
of fallout radiation on the roof of a structure. In order to do this, 
experimental work on the subject was reviewed and compared with theory. 

A number of reported experiments were examined and it was found that much 
of this work had been on structures such as houses and barracks which had 
very complicated geometrical configurations. Work reported on such struc¬ 
tures gave values of radiation exposure valid only for the structure 
tested and accurate generalization of the results to other structures 
appeared to be extremely difficult. Since the aim of this project was 
to best the method of calculation of protection against fallout radiation, 
only those experiments which had been conducted on structures of reason¬ 
ably simple geometry were used to compare with theory. The only such 
experiment reported at the time of writing was that at the Nuclear Defense 

(3) 

Laboratory . Preliminary results from work at the Protective Structures 

(4) 

Development Center were also available. 


To provide a comparison between theory and experiment, it was neces¬ 
sary to compute the exposures due to a plane isotropic source penetrating 
steel slabs, since steel slabs were used in the NDL experiment and cal¬ 
culations for steel using moments theory were not available. At the same 
time, calculations were done for concrete and for air, for a plane isotropic 
source, at several energies. 

Three factors which may cause error in comparison of theory and experi¬ 
ment were considered. These are the effect of source anisotropy in free 
field measurements, the effect of lack of source reflection in roof experi¬ 
ments, and the effect of error in estimating the thickness of the slab 
used in the experiment. The order of magnitute of errors due to the above 
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three factors were estimated and they were found generally to be small 
but not necessarily negligible, as reported in Section 3. 

(3 4) 

The experimental measurements available ’ were not corrected for 
these effects in this report. The experimental results were compared 
with the theoretical exposures and reduction factors calculated here (see 
Appendix A); agreement was often better than 10%. 

Tn the appendices, two subjects of related interest which arose 
during this investigation are reported. These are the detailed results 
of exposure distributions and reduction factor calculations in iron and 
concrete due to plane isotropic sources of various energies, calculated 
by moments theory; and a review of ground contribution to fallout radiation 
penetration through simple structures. 

2 THE LOGIC OF THE METHOD OF ANALYZING STRUCTURES FOR FALLOUT PROTECTION 

Before proceeding with a comparison of theory and experiment, it is 
helpful to know where the theory fits into engineering technology and 
how experiment relates to the theory. The relation of theory and experi¬ 
ment are discussed in this Section, together with some general considerations. 

In engineering design, the properties of a structure or other item 
are specified, and the components are then selected and proportioned so 
as to give the required results; in analysis, a structure is completely 
specified and then some property, such as its behavior under an unusual 
condition, is calculated. 

If a calculation method is available which permits direct computation 
of the size and specifications of the structure to satisfy the design 
criteria, then the structure may be directly designed. Many structures 
cannot be directly designed. The procedure then is to select a design, 
analyze it, amend the design to satisfy some given criteria, analyze again, 
and so on. We then have a feed-back interative process where analysis 
is one part of the design procedure. 

The hypothetically most accurate and direct approach to the problem 
of analyzing a structure for its ability to shield people from fallout 
radiation would be to subject it to actual fallout and measure the 
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exposure received. This Is not practical, nor is it practical to subject 
every structure to simulated fallout radiation. 

For this reason, structures are divided into classes, categorized 
by some general properties such as size, construction material, etc. 

This implies that we know which properties are significant for this prob- 
4 lem, and which are not, so that the process of generalization has small 

effect on the accuracy of the analysis. 

It is not practical to do experiments which cover all relevant regions 
of all relevant variables even within one class of structure. Recourse 
is then made to theory, which is a generalization of the practical problem 
in a way that can be applied with reasonable effort and accuracy to any 
structure in the class to which it pertains. 

Before the theory of radiation protection can justifiably be used, 
it must be established for which situations it is applicable and how 
accurate it is. This is done by comparing theory with experiment, and 
, the experiment must be as careful and as accurate as possible so that it 
gives some evidence as to the validity of the theory. Once that is 
established, the theory may be used with an established degree of confi¬ 
dence, for any particular problem. 

Since the biological effects of radiation dose are relatively insen¬ 
sitive to small changes of dose, great accuracy in the knowledge of the 
shielding capability of any one structure is superfluous. If, therefore, 
an experiment were performed on a structure simply to obtain the protec¬ 
tion factor of that structure, it need not be extremely accurate. But 
if it is performed to check a theory, it must be well understood and 
carefully carried out, so as to give limits of confidence to the theory. 

By this means a small number of careful experiments, linked to a theory, 
can replace much random empiricism. 

In order to draw conclusions about the method of calculation for 
analyzing structures against fallout radiation, it is convenient to con¬ 
sider both the experimental and practical situations, and the theoretical 
model, in five parts (see Table 1): 
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(lit The spatial distribution of the source; 

(2) The uuciear properties ol the source; 

.3) The spatial distribution of the structure; 

(4) The nuclear properties of the structure; 

(5) The definition of the dependent variable, exposure. 

These five aspects are discussed below. 

(1) It is assumed that real fallout will be uniformly distributed 

in plan. The validity of this assumption is not discussed in this report. 
The source is treated as a uniform plane in theory, and this can be 
accurately simuiaitd by experiment. 

(2) The interactions of gamma radiation with matter are well under¬ 
stood over the relevant energy range. Conclusions drawn about the theory 
on the basis of experiments with sources other than fallout may therefore 
be taken as valid for application to real situations. There is some dif¬ 
ficulty in knowing just what to consider to be the source properties in 
an experiment, on account of the fact that experimental sources are often 
large relative to the mean-free-path length of radiation within them, and 
because in calculations using moment theory to find the penetration through 
a roof slab the exposure originating from fallout on the ground (skyshine) 
is not separated from the exposure originating from fallout on the roof. 
This is discussed in Section 3. 

(3) The theory assumes the roof of the structure to be an infinite 
homogeneous slab whereas real structures consist of columns, walls, beams, 
and finite slabs. 'his is the major point of differente between theory 
and reality, and experiments should be, and are, designed to check the 
validity of the theory in different geometries of structures. 

The material geometries of the experiment correspond to those of 
real structures. In order to use experimental evidence to draw conclu¬ 
sions about the relevance of theory to reality, differences in source 
spectra between experiment and theory must be accounted for. 

(4) Nuclear properties of structural materials with respect to ganuna 
radiation are known. The effect of the small Inaccuracies in available 
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nuclear data is much less than the effect due to the relatively large 
dimensional tolerances allowed in the sizes of structural members. 

(5) The detectors used in experiment measure exposure. This is the 
dependent variable in the theory. It is worth pointing out that the 
absolute calibration error of even the best detector is of the order of 
2% or more, and some detectors used have inaccuracies an order of magnitude 
greater. 

Table 1 shows in summary form the areas of similarity and dissimi¬ 
larity as discussed above. 


TABLE 1 

Summary of relation of theory and experiment to assumed real situation. 
0 represents similarity to assumed real situation 
X represents difference compared with assumed real situation 


Property 

Theory 

Experiment 

(1) Source, Spatial 

0 

\ 

X 

(2) Source, Nuclear 

0 

X 

(3) Structure, Spatial 

X 

0 

(4) Structure, Nuclear 

0 

0 

(5) Detector response 

0 

0 


3 ANALYSIS OF SOME POSSIBLE CAUSES OF EXPERIMENTAL ERROR 
3.1 Sources and Source Calibration 

The reduction factor of a structure is expressed as the ratio of the 
exposure observed at a given point in a structure to the exposure which 
would be observed three feet above a smooth, infinite, unobstructed, 
uniformly contaminated plane. Experimentally these exposures are derived 
from separate measurements. Each exposure is normalized to an arbitrary 
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source strength, usually 1 curie/ft . and the ratio then taken. The 
protection factor is defined as the reciprocal of the reduction factor. 

If the same source were used for the structure and the f .-ee field 
expe'iment, and if it were isotropic there would be no need to normalize 
the dose to a source strength, because the exposure observed in the struc¬ 
ture divided by the free field exposure would be the protection factor. 

If an anisotropic source were used in either exp 2 riment, normalization 
of the observed exposure to a calibrated source strength would require 
that the source anisotropy be allowed for. 

In practice, the sizes of the sources and capsules used are not 
small compared with the mean-free-path length of radiation in them, nor 
are they spherically symmetrical. This leads to two effects: 

(a) Degradation of the source spectrum: 

(b> Anisotropy of the emitted exposure with respect to the source 
orientation. 

When the source calibrated in air, the low energy components of 
the spectrum are detected. When the source is placed on a barrier, the 
low energy components are more rapidly attenuated than the higher energy 
components, If the results are compared with a theory derived for the 
source energy, the theory will overestimate the transmitted exposure. 

(5 

Berger and Doggett compared thetr calculations with an experi- 

(61 

ment of Kirn, Kennedy and Wyckoff and corcluded that neglect of the 
energy degradation from a 1 cm thick Co-60 source gave rise to an 
ovemalculation in exposuie of 18 . for barriers approximately 

120 psf thick. 

The error introduced in an experimen* by neglecting the degraded 

source spectrum can be estimated ll the source spectrum is known, by 

using the penetration calculations for various incident energies given 

in this report in Appendix A. Unfortunately the source spectra were 

not determined lor sources used in the experiments. This effect may 

tn fact turn out to be small in some experiments but no firm evidence 

« 

is available to allow this assumption. 

* 

A program to calculate by Monte Carlo methods the spectr as of gamma 
radiation emitted from a source is at the present time being written, 
and when available l* should help to solve this problem. 
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The experimental value of exposure in air three feet above the source 
plane is obtained by measuring the exposure, either when the source is 
pumped through a tube on the ground, or by placing a detector in many 
locations when a source is placed in one location and the results summed 
to Simula* a plcne field. An estimate of the effect oi source anisotropy 
in the latter procedure is made in Section 3.2. 

3.2 Effect of Source Anisotropy in Free Field Measurements 

(7) 

Th' free field measurements made at the Nuclear Defense Laboratory 
had been carried out by placing detectors in a number of positions relative 
to a single source on a field and summing the observed exposures to s-nru- 
latt: an infinite plane source. Since the source was probably not iso¬ 
tropic, the observed exposure, when normalized to the calibrated source 
strength, may have been in error. An estimate of the error due to this 
effect was made as follows. 

Since the anisotropy of the source radiation was not known for the 
- ources used at the Nuclear Defense Laboratory, an estimate of the order 
of magnitude of the effect was made for a hypothetical case, using the 
reported exposure distribution around a Co-60 source, given by Burson 
and Summers' and shown in Fig. 1. 

The reduction factor due to area n of Co-60 in a field is given by 


~ - P(d) A/4rrd 2 3.2.1 

o 


where A and d are shown in Fig. 2 and P(d) for Co-60 is given in 
Fig. B-19 of the monograph by Spencer Referring to Fig. 2, this can 

be written for a plane isotropic source 


D 

D 



P(d) 2nr dr/4nd“ 


3.2.2 


The exposure emitted by the Co-60 source shown in Fig. 1 can be written, 
with some approximation as 
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J i 1 - 0.0033 G ) . 

[ °- 5 
Including this anisotropy effect in eq. 


-60° 

< 9^60° 

3.2.3a 

60° 

< 9^90° 

3.2.3b 


3.2.2 gives an approximate relation 



D 

D 





F(6> P(d) r dr/2d 2 


3.2.4 


where D is now the exposure at the detector due to a uniform infinite 
contamination corresponding to a calibrated source intensity ; the 
source energy is here assumed to be 1.25 MfeV with no degradation. 

Changing varietbles, ttie relative exposure due to an infinite plane 
source is 


D 

D 




Fee' 


P(d) 


d(sin9> 

sin9 


3.2.5 


To find the effect of source anisotropy, the ratio 





F(e: p d'< 


d t sin9) 
sin6 


P(d> 


d i sj n 9, 
sin9 


was calculated and found to be 0,96 . 

This calculation shows that if iree field measurements were acne by 

(7) 

the method described by Fexroad and Schmoke using a source similar to 

<81 

that described by Burson and Summers , source anisotropy would decrease 
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the measured relative exposure by approximately 4% . If the source were 
such that the exposure on the axis were greater than that observed normal 
to the axis, source anisotropy would increase the relative exposure. 

3,3 Effect of Lack of Source Reflection in Roof Experiments 

The barrier and geometry factors in present methods of calculation 
of fallout protection are calculated usually by the method of moments. 
These assume, among other things, an infinite extent of material on that 
side of the source which is away from the detector. 

In reality, fallout would be in a thin layer on and around a struc¬ 
ture with air above it. If the fallout field were infinite in extent, 
or greater than three mean-free-path lengths (about 1500 feet) in radius, 
and if the building height were small compared with a mean-free-path 
length in air, then the contribution of radiation scattered from the air 
above would equal the theoretical contribution due to the radiation 
scattered back by the hypothetical infinite medium. 

If the source in a roof experiment were covered with an infinite 
thickness of the same material as the roof, the experimental configura¬ 
tion would be identical with the theoretical source configuration in the 
moments calculations,and this would include the equivalent of "skyshine." 
Since this was not always the case in past experiments, an estimate must 
be made of how much the theory may have differed from the experimental 
results on account of the experimental source characteristics. 

(5) 

Berger and Doggett calculated buxld _ up factors for finite and 
semi-infinite barriers for plane normal sources. The difference between 
the results for the two barriers, for a 1 MeV source on iron, at one 
mean-free-path, was 1% . For an isotropic source, the difference may 
be expected to be larger. 

An approximate calculation can be made as follows: From the results 
presented in Appendix A the exposure penetrating into iron from a i.25 
MeV plane isotropic source is known. We wish to subtract the exposure 
due to the radiation which moves initially away from the detector and 
then cresses the source plane to be detected. 
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SUMMARY 


A comparison is made between experiment and theory for the calcu¬ 
lation of shielding by a structure against fallout radiation. 

Experimental results for exposure penetrating a roof slab, and 
for reduction factor, are found to agree with moments theory calcu¬ 
lations, often to better than 10% when the geometry factor L (u),X) 
was used. 

A comparison between experiment and theory may be inaccurate due 
to anisotropy of the experimental source, to lack of source reflection 
in root experiments, and due to error in estimating the thickness of 
the roof. The magnitudes of the errors due to these effects are inves¬ 
tigated, and found to be small, but not necessarily negligible. 

Detailed results of penetration in iron and concrete, due to plane 
Isotropic sources of various energies are given, and a review of ground 
contribution to fallout radiation penetration through simple structures 
is appended. 
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Berger and Raso give the following values for the number albedo 
A from a 1.25 MeV isotropic source: on iron 0.295 ; on concrete 
0.383. 

(9) 

Berger and Raso calculated the reflected energy spectrum ior a 
1 MeV isotropic source on iron. This is shown in Fig. 3. On the basis 
of this result, i‘ will be assumed that all the reflected photons due to 
a Co-60 source on iron or concrete have energy 0.25 MeV . 

The current due to scattered radiation which crosses the source 
plane will be A times the source current and is assumed to be at 0.25 

MeV . Assuming, without justification, that this is isotropic at the 
source plane, the penetration of the scattered radiation can be found 
from the moments calculation at 0.25 MeV . The exposure at depth Z 
psf in the barrier, due to reflection above the source, will be 

A • D(Z, 0.25 MeV) 

D(Z, 1.25 MeV) 

This calculation will be approximate, because the reflected photon dis¬ 
tribution will be more peaked than isotropic and because the energies 
of the reflected photons are not at a single energy, but an estimate of 
the order of magnitude of the effect is obtained from this calculation. 
The calculated ratios of exposure reflected by source cover to exposure 
in an infinite medium, are shown in Table 2. 

The above calculation gives approximately the same result as cal¬ 
culations using the functions given by Spencer ^\ In his notation: 

S<"d) is the ratio of exposure due to radiation from the upper hemisphere 
( akyshine") to the reference exposure 3 feet in air above the source 
plane; 

S'(X) is the ratio of exposure at X psf to the exposure at 0 psf, due to 
a plane source in an infinite medium, emitting radiation Isotropically 
into the half-space in which the detector is buried; 

L(X) is the rsito of exposure due to a plane infinite isotropic source, 
at X psf in an infinite medium, to the reference exposure 3 feet in air 
above the source plane. 
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A comparison between the two methods of calculation is obtained by 

. A • D(Z, 0.25 MeV) , _ , , , _ 

comparing the value of - —: ~rr - in the previous calculation 

r/h's R« ✓y'* D ^ z > 1,25 MeV ' 


with 


S(d) S'( X) 
L(X) 


in this calculation, when d = 0 . Note that both 


calculations refer to infinite source planes; solid angle fraction uo 


is then one and S (X,uj) = 1 . 

a 

The latter approach was used by Raso and Woolf 


(17) 


when comparing 


their Monte Carlo calculations of roof contributions with the moments 
calculations of Spencer^. When the quantity S(d) S'(X) S (<ju,d) was 

H 

subtracted from the value of L(X) L (X,u)) , good agreement was found 

& 

with the Monte Carlo results. 


TABLE 2 


Ratio of exposure reflected by source cover to total exposure for iron 
and concrete. Results from approximate calculation, using albedo and 
penetration curves; and using curves of NBS 42^'. 


Material 

Z PSF 

Calculating 

Procedure 


25 

50 

100 

150 


Iron 

0.025 

0.011 

0.003 

0.0005 

mm 






D(Z, 1.25 MeV) 

Concrete 

0.062 

0.037 

0.013 

0.004 


Concrete 

0.077 

■ 

0.031 

— 

S(d) S*(X) 

L(X) : d=0 


The results of the two mett of calculstion are given for concrete 
in Table 2 and show reasonable agreement at low barrier thicknesses, when 
the correction is more significant; at greater thicknesses agreement is 
not as good, but the correction is then not important. 

The calculations show that lack of source reflection may account for 
about 3% of the exposure in iron and 6% in concrete at 25 psf, about 1% 
in iron and 4% in concrete at 50 psf. and less at greater distances. This 
difference between theory and experiment may be avoided simply by having a 
reflector over the source. 
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3.4 Errors in Estimating Roof Thickness 

During the course of this investigation, a question arose as to the 
error provoked when relating theoretical to experimental exposures, by 
assuming that a slab has exactly the nominal thickness when in fact it may 
not have. An extimate of this cause of error is made as follows. 

Consider a roof slab, upon which is placed a uniformly distributed 
plane isotropic source. The radiation exposure within the slab can be 
approximately given by D = D^e ^ when t is the slab thickness and 
li is the effective broad beam attenuation coefficient for the slab. The 
relation between errors in D and t is given by —^ fit . 3.4.1 

The effective broad beam attenuation coefficient for a plane iso¬ 
tropic Co-60 source on steel is found from Fig. A2 of Appendix A to be 
0.027 psf 1 , or 1.1 in. 1 Assuming this is known with zero error 
—= -1.1 6t where t is in inches. 3.4.2 

The difference between real and nominal thicknesses is demonstrated 
by observations made during an unpublished experiment at the University 
of Illinois. Sheets of glass, graphite, aluminum, and lucite, each 

4.5 sq. ft. in area and nominally l" thick, were measured on a 6" sq. 
grid. The mean dimensions found are shown in Table 3. 

TABLE 3 

Measured thickness of nominal i'' muterial. The mean of 18 micro¬ 
meter readings is given 


Materia] 

Thickness Cinches) 

A1uminum 

0.980 

Class 

1.005 

Luci*e 

0.967 

Graphite 

1.002 


The American and Iron Steel l 
for structural plates to be minus 


( 12 ‘* 

nstitute gives the rolling tolerances 
0% plus 7% . If we assume that steel 
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used in an experiment Is 
then 6t =0.03 x t and 

-j p ~ = 1.1 x .03 x t 
= 0,03 x t 


3% thicker than the assumed nominal value, 
the observed exposure will be decreased by 


3.4.3 


The observed exposure will then decrease by about 3% for each inch 
thickness of steel. 


4 COMPARISON OF THEORY AND EXPERIMENT 

(3) 

4.1 The Nuclear Defense Laboratory Blockhouse Experiment 

In this experiment sources of Co-60 and Cs-137 were placed on 
the roof of a blockhouse so as to simulate a plane source; detectors were 
placed within the blockhouse and exposures ure measured. In a separate 
experiment, the same sources were placed » en field and the exposure 

three feet above the effectively infinite , source measured. Dividing 
the former exposure by the latter gave the reduction factor. 

The blockhouse had internal dimensions 12 ft by 12 ft by 8 ft high. 
The walls were of concrete, 4 inches thick. The roof, which was supported 
by the walls and by a joist at the blockhouse center line, consisted of 
1/2 inch plywood, and various thicknesses cl steel or concrete above this. 
The source rested on a 1/2 inch thickness of lucite on the roof. 

The detectors were Victoreen Model 208 lOmr dosimeters. These were 
calibrated to ±3% against a Victoreen Model 130 r-meter which itself 
had been calibrated to ±2% at the National Bureau of Standards. Twenty- 
four detectors were placed within *he blockhouse at various positions. 

The sources used were 0.34 curios of Co-€0 and 1.35 curies of 
Cs-137. The cobalt was metal in a stainless steel capsule; the cesium 
was CsCl also in a stainless steel capsule. The sources were calibrated 
by holding them a known distance from NBS-calibrated ion chambers. The 
calibration of the ion chambers was accurate to ±2% , and it is expected 
that the calibration of the source relative to the ion chambers was 
accurate to about 2% , so the source calibration was probably accurate to 


( 0 . 02) 2 + ( 0 . 02) 2 


3 % 
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Each source was placed on the roof in a number of locations and the 
exposure rates thus observed were summed so as to be equivalent to a 
uniform source distribution on the roof. 

The experimental exposures at various heights on the blockhouse 
center lines were compared with theoretical values by the NDL workers, 
in the following way: 

(a) The exposures observed with a steel roof were multiplied by a 
factor of 1.11 to make them equivalent to those which would be observed 
with a concrete roof. The factor 1.11 was used as the average ratio 

of ordinates of the curves of exposure vs effective mass in concrete 
and iron for a 1 Me» plane normal source (Fig. 22.1 of reference 1). 

(b) The reduction factor was calculated by dividing this exposure 

rate by the exposure rate D , observed three feet above an infinite 

o 

contaminated field. The values of D were observed in previous work 

(7) ° 

at the same site. 

(c) The theoretical reduction factor was taken from Fig. 28.19 of 

NBS 42^\ which gives values of L(X) L (X,uj) calculated for water and 

a 

a 1.12 hr fission source. 

Comparison of those experimental results with theory was repeated 
in this investigation, making use ot the newly compu’ed penetrations in 
steel (see Appendix A' and using a procedure as follows: 

(a) Roof thicknesses were expressed as actual pounds per square foot 
Z. not as equivalent mass thickness X isee Appendix A). 

tb' The exper.mental exposure rates for steel roofs given in Tableo 

2 

3.1 and 3.2 of the '*DL report in roentgen/hr per curie/ft were divided 
by 1.11, so as to eliminate the correction’' previously made. The sources 
had been calibrated using specific gamma ray constants of 14.3 roentgen/hr 
at one foot from a curie of Co-63 and 4.2 roentgen/hr at one foot 
from a curie of Cs-13 7, '’he values corresponding to 86.9 roentgen 
per erg/gw are 14.0 roenigen/hr at one foot from a curie of Co-60 and 
3.4 roentgen/hr at one foot from Cs-137. The values of exposure rate 
per curie, used in the theory and injthe reduction of experimental data, 
were made equal. 
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(c) The exposures obtained from the moments calculation for 1.25 

2 

MeV were multiplied by 5284 to convert from MeV/gm per photon/cm to 

2 

roentgen/hr per curie/ft . This constant is obtained as follows: 

1 MeV = 1.602 x 10" 6 erg 
1 roentgen = 86.9 erg/gm 

1 curie of cobalt = 7.4 x 10^ photons/sec 

2 2 
1 ft = 929.03 cm 

2 

Thus, 1 MeV/gm sec per photon/cra sec 


1.602 x 1 0 
86.9 


-6 


1U 


7 „4 x 10 
* 929.03 

= 5284 roentgen/hr per curie/ft' 


x 3600 


Since cesium-137 emits 0.84 photons per disintegration whereas 

cobalt-60 emits two photons per disintegration, the conversion factor for 

2 

cesium-137 will be 2219 roentgen/hr per curie/ft . 

(d) The geometry factors L (qa,X) and L (uj,X) were obtained from 

a (1) C 

Figs. B27, B28, B29, and B30 of NBS 42 , which are for Co-60 and Cs-137 

in concrete. Geometry factors have rot been calculated for steel, because 
it is expected that these factors will be much less dependent on the 
difference between mass attenuation coefficients for steel and concrete 
than the barrier factor. 


(e) The theoretical exposure rates were computed from D(Z)L <a),X'' 

a 

and D(Z)L (tu,X) , where D(Z' was obtained as shown in (c) above, 
c 

These are compared with the NDl observations in Table 4. 

(f) The reduction factor was ca.culated by dividing the exposure 
rate obtained in te) above by the calculated exposure in the air three 
feet above the ground as shown in Appendix A. These calculations are 
compared with the NDL olservations in T able 5. 


Tables 4 and 5 and Figures 4 though 7 show the calculated and observed 
exposure rates and reduction factors. The curves are not plotted beyond 
45 paf because the observations at 50 psi were for a concrete roof, whereas 
the other observations and calculations are for ste?l. 

There are several possible causes of error in comparing theory to 

(3) 

experiment; since these were not all evaulated in the report an order 
of magnitude estimate of these errors is made here. 
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TABLE 4 


Ratio of Theoretical Exposure Rate 
2 

(mr/hr per curie/ft ) to Observed Exposure Rate 


Source 

Detector 

(11 

Barrier Thickness PSF 

Theoretical 
Geometry 
Factor Used 

Location 


18.0 + 

33.7* 

44. 9 + 

EES■ 

Co-60 

C 6 

0.73 

1.08 

1.08 

1.09 

0.98 

L (uu,X) 
a 


C 3 

0.41 

0.97 

1.03 

1.08 

0.97 



C 0 

0.24 

0.74 

0.83 

0.92 

0.83 



C-l 

0. 20 

0.79 

0.86 

0.85 

0.78 


Co-60 

C 6 

0.73 

1.17 

1.26 

1.29 

1.22 

L c (-jj,X) 


C 3 

0.41 

0.95 

1.07 

1.18 

1.17 



C 0 

0.24 

0.72 

0.83 

0.JH 

0.94 



C-l 

0.20 

0. 7S 

0.84 

0.91 

0.90 

i 

Cs-137 

C 6 

0.73 

1.06 

1.12 

1.26 

1. 13 

L (x. X) 
a 


C 3 

0.41 

1.06 

1.02 

1.06 

1.03 



C 0 

0.24 

0.97 

0.91 

0.92 

0.90 



C-l 

0.20 

0.98 

0.87 

0.91 

0.94 


Cs-137 

C 6 

0.73 

1.26 

1.37 

1.55 

1.45 

L (x.Xl 
c 


C 3 

0.41 

0.96 

1.25 

1.35 

1.43 



C 0 

0.25 

0.81 

1.03 

1.11 

1.22 


■ —. 

C-l 

0.20 

0.78 

0.97 

1.06 

1.18 



* Steel was used in the experiment for these thicknesses. 

* Concrete wss used in the experiment for this thickness only. 
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TABLE 5 


Ratio of Theoretical Reduction Factor 
to Observed Reduction Factor 


Source 

Detector 

Location 

(D 

Barrier Thickness PSF 

Theoretical 

18. 0 + 

33. 7 + 

44. 9 + 

50.3* 

Geometry 
Factor Used 

Co-60 

C 6 

0.73 

1.00 

1.00 

1.00 

1.02 

L (u),X) 


C 3 

0.41 

0. 82 

0.95 

1. 00 

1.03 



C 0 

0.24 

0.69 

0.77 

0.84 

0.85 



C-l 

0.20 

0.71 

0.74 

0.-81 

0.78 


C0-60 

C 6 

0.73 

1.09 

1.17 

1.20 

1.26 

L (u),X) 








c 


C 3 

0.41 

0.80 

0.98 

1.09 

1. 22 



C C 

0.24 

0.67 

0.77 

0.91 

0.96 



C-l 

0.20 

0.71 

0. 74 

0.85 

0.91 


Cs-137 

C 6 

0.73 

0.94 

1.00 

1. 13 

1.10 

L a (u),X) 


C 3 

0.41 

0.94 

0.91 

0.97 

0.99 



C 0 

0.24 

0.87 

0.81 

0.82 

0.87 



C-l 

0.20 

0. 86 

0.79 

0.84 

0.94 


Cs-137 

C 6 

0.73 

1.12 

1. 2? 

1.39 

1.42 

L c (tu,X) 


C 3 

0.41 

0.86 

1.11 

1.20 

1.42 



C 0 

0. 24 

0.71 

0.91 

1.00 

1.17 



C-l 

o.?o 

0. 70 

_ 

0.86 

0.95 

1.17 



Steel was used in the experiment for these thicknesses, 

* 

Concrete was used in the experiment for this thickness only. 
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(a) The error ia the calculated values was of the order of 2% due 

to difficulties in reading accurately the values of L (u),X) and L (d,X) 

( 1 ' 3 C 
from the graphs in N?S 42 '; 

(b) Uncertainty in absolute calibration of the source strength is 
expected to have been approximately 37. (see text); 

(7) 

vC) Because of anisotropy of the source in free field measurements , 

2 

D for Co-60, when expressed as exposure rate per curie/ft , may possibly 
o 

have been in error by approximately 4% (see Section 3.2); 

(d) If the Co-60 source were entirely unreflected then the theoretical 
exposure rate would have been higher than the experimental exposure rate 

by approximately 2% at 25 psf and 1% at 50 psf in steel (see Section 3.3). 
Since in fact the finite size of the source would have caused some radia¬ 
tion to be reflected back into the roof this effect may have been small, 
for the steel roofs tested. For concrete, the effect would have been 
larger, but still small; 

(e) The nominal thicknesses of the steel roofs appear to have been 
used for the calculations in the NDL report. II is possible that the 
thicknesses were in fact greater than nominal. If they were 3% thicker 
the exposures would have been cvercalculated by approximately 1% at 18 psf, 
2% at 34 psf and 3% at 45 psf for the Co-60 source (see Section 3.4); 

(f) Because of degradation of the source spectrum, the measured 

2 

values of exposure rate per curie/ft through the barrier would have been 
lower than the theoretical values. ,T, he magnitude of this effect is, for 
the present, unknown (see footnote to page 6'. 

Ihe effect listed in f cj above would have reduced or increased the 
free field exposure rate, depending upon the geometry of the encapsulated 
source, whereas those listed in (d), (e), and (f) above would have decreased 
the exposure rate observed in the blockhouse. This discussion of possible 
errors is based on general considerations and no judgement can be given 
here as to whether ~hese causes of difference between theory and experi¬ 
ment were actually present. 
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Agreement between calculated and observed reduction factors appears 

to be better for larger solid angles and greater barrier thickness, and 

better agreement is found between theory and experiment when the geometry 

factor L (cd,X) is used than when the factor L (oj,X) is used, 
a c 

The geometry factor L (u),X) is calculated for an infinite plane 

source in an infinite medium and a limited solid angle of detector response 

the factor L (u),X) is calculated for a circular plane source in an 
c 

infinite medium and a detector response over 4tt steradiams. Since in 
fact the barrier has a finite thickness, and the detector is some dis¬ 
tance from the barrier with only air in its close vicinity, the factor 
L (u>,X) may indeed be expected to give more accurate results. 

Of the 16 values of reduction factor calculated for Co-60 (using 

L (u),X) ): 

8 

7 fall within ± 10% of experiment; 

11 fall within ± 20% of experiment; 

15 fall within ± 30% of experiment. 

For Cs-137 (using L (u),X) ): 

8 

8 fall within ± 10% of experiment; 

15 fall within ± 20% of experiment. 

( 2 )* 

4.2 The Protective Structures Development Center Idealized Roof 

In this experiment, a circulating source was placed on the ground 
with concrete slabs above. Detectors were placed at various heights above 
the slabs. 

The reduction factors were calculated by dividing the observed expo¬ 
sure rates by the exposure rates previously measured three feet above 
an infinite plane of contamination. These experimental reduction factors 
were compared with the theoretical values of L(X)L (X,ua) and L(X)L (X,u>) 

A C 

and this comparison is repeated here by taking the values given in that 

♦ (4) 
report. 


* 


At the time of this writing, preliminary results only are available. 







20 


The L(X)L (X,uj) curves were closer to the observed results than 
a 

the L(X)L (X,uj) curves, at all values of X . The differences are 
c 

listed in Table 6, 


TABLE 6 


Comparison of Theory and Experiment for the PSDC Idealized Roof. 
Ratio of Calculated Reduction Factor to Observed Reduction Factor 


Barrier 

PSF 

U) 

Theoretical Value 

L(X)L Cut) i X) 
a 

L(X)L (uu.X) 

C 

0 | 

0.02 

■. 

1.00 

0.95 


0. 10 

1.00 

0.94 


0.60 

1.00 

0.94 

48 

0.02 

1.12 

1.19 


0.10 

1.15 

1.25 


0.60 

1.17 

1.43 

96 

0. 02 

1.07 

1.40 


0. 10 

1.13 

1.47 


0.60 

1.09 

1.36 

144 

0.02 

1.00 

1.50 


0. 10 

1.00 

1.55 


0. 60 

1.00 

1.35 


In this experiment, the unprotected exposure was observed with detec¬ 
tors above the circulating source. The protected exposure was obtained 
by inserting concrete slabs between source and detectors. In this 
arrangement, error duo to lack ol source reflection was eliminated since 
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the ground served as a reflector. The source was a cylinder which cir¬ 
culated so that the circular suiface faced the dosimeters for a large 
portion of the travel. Since source anisotropy is usually most marked 
beyond the ends of the source, it is expected that this arrangement 
gave less error due to source anisotropy than the previously discussed 
experiment where a single source was placed on the ground or on the roof. 

Because of the two above-mentioned reasons, and because the experi¬ 
ments w. re compared with calculations for Co-60 incident on concrete 
slabs, th^-.e experimental results were compared with existing theo¬ 
retical calculations, no further calculations being carried out. 

5 CONCLUSIONS 

Various possible causes cf error in comparing experiment with theory 
were onsidered in this investigation. The conclusions may be summarized 
as follows: 

(1) Source anisotropy in free field measurements using a stationary 
source may cause an error on the order of 4% when normalizing the 
measured exposure to the calibrated source intensity. It must be empha¬ 
sized that this is an approximate result only. If more accuracy in 
estimating this cause of error is required, more detailed calculations 
must be made for any particular experiment. 

(2) The lack of reflector above a source in a roof experiment will 
cause larger error for a concrete than for a steel roof; this error for 
concrete will be of the order of 6% at 25 psf, 3% at 50 psf, and 
decreases to less than 1% above 100 psf. The source should have a 
reflector so as to include the equivslent of "skyshine" which is implicit 
in the moments theory calculations. 

(3) Errors in the assumed value of thickness of a steel roof will 
cause an error in the calculated exposure of about 1% per inch thick¬ 
ness, for an error of 1% in assumed roof thickness. 

(4) Comparison of the results of experiments with theoretic; 1 values 
often show agreement to within 10% when the geometry factor L^fui.X) is 
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used. There is a trend for the theory to predict lower values of expo¬ 
sure and reduction factor as solid angle decreases. 

(5) It is recommended that the geometry factor L a (uu,X) be used 
for roof penetration calculations in simple structures, and in that case 
an accuracy of greater than 10% in predicted reduction factor may often 
be expected. 
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Enlarged view 
of source 


SCHEMATIC VIEW OF FREE FIELD EXPERIMENT AT 
THE NUCLEAR DEFENSE LABORATORY (REFERENCE 7) 


FIG. 2 


2 


FIG. 3 ENERGY SPECTRUM OF I MeV PHOTONS 
REFLECTED FROM IRON 


ISOTROPIC INCIDENCE, I MeV SOURCE 
FROM TABLE A 19, BERGER S RASO 

NBS 5982, I960 







FIG.4 COMPARISON OF EXPERIMENTAL 
AND THEORETICAL EXPOSURE 
RATES 


EXPERIMENTAL — SCHMOKE AND REXROAD 1962 

o Detector Position C6 

• C3 

a CO 

• C-l 

Theoretical - L a (w,X)D(X) 

- L c (w,X) D (X) 


Source Co-60 
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FIG.5 COMPARISON OF EXPERIMENTAL 
AND THEORETICAL EXPOSURE 
RATES 


EXPERIMENTAL - SCHMOKE AND REXROAD 1962 

o Detector Position 06 

• C3 

* CO 

» C-l 

L„ (asX)D(X) 

L c (u>,X) D (X) 

Source Cs-137 


Theoretical 



o 


C6 









REDUCTION FACTOR 


30 



BARRIER PS F 



FIG.7 COMPARISON OF EXPERIMENTAL 
AND THEORETICAL REDUCTION 
FACTORS 


EXPERIMENTAL - SCHMOKE AND REXROAD 1962 

o Detector Position C6 
• C3 

a CO 

■ C-l 

Theoretical -- L a (<u,X) D (X)/D 0 

—-L c («,X) D (X) /D 0 

Source Cs-137 


X 
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APPENDIX A 

CALCULATIONS OF EXPOSURE AND REDT.'CTION FACTOR DUE TO 
PLANE ISOTROPIC SOURCES IN IRON AND CONCRETE 


Exposure distributions were calculated in infinite media of concrete, 
iron and air for plane isotropic sources, by the method of moments. Use 
was made of the code available at the National Bureau of Standards for 
this. Source energies were 1.33, 1.25, 1.17, 0.6616, 0.25, 0.10, and 
0.05 MeV. The material cross sections were taken from those published 
by Grodstein^*^ and McGinnies^^. 

Table A1 gives the exposure in iron and in concrete at a number of 
depths for the various source energies. These values are plotted in 
Figures A1 to A4. 


Table A2 shows the following quantities for each of five source 

energies; the exposure in air three feet from the source plane in MeV/gm 
2 

per photon/cm , and the ratio of scattered to total exposure at the 
detector. 

2 2 

Conversion from MeV/gra per photon/cm to roentgen/hr per curie/ft 
is shown in Section 4.1. Reduction factors may be obtained by dividing 
the exposures in concrete or iron by the exposure in air three feet 
from the source plane at the relevant source energy. 


It should be noted that the exposures are given as functions of the 

mass thickness Z . This is not to be confused with the equivalent 

(1 ^ 

mass thickness X of NBS 42 





33 


TABLE A1 


Exposures in MeV/gm per photon/cm for Iron and Concrete 
at Different Source Energies and Depths 


E 

o 

MeV 

p z 
o 

Iron 

Concrete 

z 

psf 

Exposure 

z 

psf 

Exposure 

1.33 

0.5 

19.9 

0.0185 

18.6 

0.0199 


1.0 

39.7 

0.00932 

37.1 

0.0103 


1.5 

59.6 

0.00520 

55. 7 

0.00588 


2.0 

79.5 

0.00302 

74.2 

0.00347 


2.5 

99.3 

0.00179 

92.8 

0.00209 


3.0 

119 

0.00107 

111 

0.00127 


4„ 0 

159 

0.000395 

148 

0.000473 


8.0 

318 

0.00000779 

297 

0.00000769 

1.25 

0.5 

19.3 

0.0178 

18.0 

0.0193 


1.0 

38.5 

0.00901 

35,9 

0.0101 


1.5 

57.8 

0.00504 

53.9 

0.00575 


2.0 

77.0 

0.00293 

71.9 

0 00340 


2.5 

96.3 

0.00174 

89.4 

0.00205 


3.0 

116 

0.00105 

108 

0.00125 


4.0 

154 

0.000387 

144 

0.000469 


8.0 

308 

0.00000771 

288 

0.00000973 

1.17 

0.5 

18.6 

0.0171 

17.4 

0.0186 


1.0 

37.2 

0.00867 

34.8 

0.00973 



55. 9 

0.00486 

52.2 

0.00559 



74.5 

0.00284 

69.5 

0.00332 



93.1 

0.00169 

86.9 

0.00201 


3.0 

112 

0.00102 

104 

0.00123 


4.0 

149 

0.000376 

139 

0.000463 


8.0 

298 

0.00000759 

278 

0.00000977 

.6616 

0.5 

14. 1 

0.0113 

13.3 

0.0130 


1.0 

28.2 

0.00583 

26.5 

0.00709 


1.5 

42. 2 

0.00332 

39.8 

0.00421 



56.3 

0.00197 

53.0 

0.00258 



70.4 

0.00119 

66.3 

0.00160 



84.4 

0.000723 

79.5 

0.00100 


4.0 

113 

0.000274 

106 

0.000395 


8.0 

225 

0.00000597 

212 

0.00000961 
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TABLE A1 (Continued) 


0.25 


0.10 


0.05 


Iron 

Concrete 

z 

psf 

Exposure 

...... . .. 

z 

psf 

Exposure 

8.66 

0. 00400 

8.93 

0.00584 

17.3 

0. 00200 

17.9 

0.00335 

26.0 

0.00111 

26.8 

0.00206 

34. 7 

0.000636 

35.7 

0.00130 

43.3 

0.000374 

44.6 

0.000822 

52.0 

0.000223 

53.6 

0.000525 

69.3 

0.0000806 

71.4 

0.000215 

139 

0. 00000155 

143 

0.00000588 

3.00 

0.000901 

6.03 

0.00193 

5.99 

0.000389 

12.1 

0.00105 

8.99 

0.000189 

18.1 

0.000613 

12.0 

0. 0000978 

24.1 

0.000366 

15.0 

0. 0000522 

30.1 

0.000221 

18.0 

0. 0000286 

36.2 

0.000135 

24.0 

0.00000893 

48.2 

0.0000507 

47.9 

0.000000110 

96.4 

0.00000108 

0.569 

0.000577 

2.92 

0.000972 

l.H 

0.000232 

5.84 

0.000448 

1.71 

0.000107 

8.76 

0.000228 

2.28 

0.0000530 

11.7 

0.000121 

2.85 

0. 0000273 

14.6 

0.0000665 

3.42 

0.0000144 

17.5 

0.0000371 

4.55 

0.00000425 

23.4 

0.0000120 

9.11 

0.0000000443 

46.7 

0.000000163 














35 


TABLE A2 

Exposures in Air Three Feet from a Plane Isotropic Source 


E 

o 

MeV 

Exposure 

2 

MeV/gm per photon/cra 

_ ^ Scattered Exposure 

Ratio —--- 

Total Exposure 

1.33 

0.0957 

0.172 

1.25 

0.0911 

0.177 

1.17 

0.0865 

0.181 

0.6616 

0.0530 

0.231 

0.05 

0.00624 

0.503 


From the above table, it can be calculated that: 

For Co-60, exposure rate three feet above source plane a 482 roentgen/hr 
2 

per curie/ft ; 

For Cs-137, exposure rate three feet above source plane = 118 roentgen/hr 
per curie/ft^ . 


The assumed density of air 


0.001205 g/cc . 
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APPENDIX B 

A REVIEW OF GROUND CONTRIBUTION TO FALLOUT RADIATION 
PENETRATION THROUGH SIMPLE STRUCTURES 


B. I Introduction 


There has been a flurry of experimental efforts in the past few 
years to check the validity of the functions used in computing the ground 
contamination contribution to fallout radiation penetration into simple 
structures, under infinite field conditions. In many cases, it is dif¬ 
ficult to design experiments which will isolate the function to be 
studied; and therefore the results obtained do not always lead to unam¬ 
biguous conclusions. It is believed appropriate at this time to review 
some of the basic technology and the conceptual basis therefor. This 
will provide an insight into the validity of previous experimental pro¬ 
grams and will be a basis for recommendations as to the directions for 
fruitful research in the near future. 


Let us consider the basic formula presently used to calculate the 
ground contribution to the attenuation function of a simple, rectangular- 
one-story structure under inf in;:v fallout field condi*: ions: 


R r 


B 

w 


C d E d ° 


S > 

w 



s . 


( 1 ) 


The notation used is roughly similar to that in the OCD Engineering Manual 

* 

and should be readily identified bv anyone familiar with the Manual. 

We will attempt »o explain this formula in ter*.* of ‘'operational'* 
concepts, that is, In terms of experteen's (perhaps Idealized) devised 
to measure them. Concepts which cannot be explained in such a manner 


This is essentially the sameas the formula lor C given in Example 
4-3 of the Engineering Manual . Our terms 0^ ind arc actually 

split into two terms in the Manual. (*ur term E rf does not appear in 
the Manual formula, but this is explained below. The remaining corres¬ 
pondence of terms is obvious. 


(13) 
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cannot be determined experimentally. The various terms in the formula 
can be related to ratios of exposure rate measurements carried out foi 
various structural configurations. These are depicted in Figure B-l 
and are discussed below. 


B.2 Explanation of Situations Depicted 

Before discussing the interrelationships of these configurations, 

or "situations," presented, it is desirable to describe them and to 

indicate the present theoretical (mathematical) approach to finding the 

exposure rates in the designated detector locations. The formulations 

( 2 ) 

are taken or inferred from Eisenhauer's NBS Monograph 76 , with further 

reference as necessary to Spencer's NBS Monograph 42^\ (These will 
hereinafter be called NBS-76 and -42, respectively.) 



in terms of the ^-function defined in NBS-42. This formula is exact, 
and in fact the ^-function is normalized so as to make equal to 
unity. 


Situation 2 provides an exposure rate behind two walls, very close 
together. The wails are of indefinite height and length, because the 
detector is so close to the walls that radiation can be assumed to be 
entering the detector from all directions (4tt steradians). The wall mass 
thickness is the same as that of the given blockhouse structure. The 
exposure rate in this situation is based on the NBS-42 function for 
penetration through a single wall (Eq. 27,7): 
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This function must be multiplied by 2, to account for contribution through 
two walls rather than one. “"he result is then slightly too high for 
small values of the mass thickness, X, because the function W always 
includes the effect of an infinitely thick backscatterer behind the detec¬ 
tor instead of a finite wall. N3S-76 indicates that this can be allowed 
for by deliberately lowering the plotted curves for small X , so that 
the resulting exposure rate is unity for X - 0 and H = 3 feet . Thus, 
D 2 ^ 2W(X,H) (4) 

for all but very small values of X . However, there are some minor 
approximations inherent in the above formula for W , so that Equation 
(4) is not rigorously correct, even for very large values of X . (See 
NBS-42, p. 3d' 

Situation 2-0 is the same as Situation 2, except that the mass thick¬ 
ness of the walls is zero. The detector response is the same as it would 

be in the open, that is, 
i 

°2-0 " I £ ,F . cose ' d(cos0) = L(K' (5) 

J-l 

where the ^.-function is defined in NBS-42 ^Equation 27.2). D 

<Ct mm U 

obviously equals when H equals three feet. 

Situation 2-°° is also similar to Situation 2, except that the mass 
thickness of the walls is approaching infinity. This is to be regarded 
as being brought about by maintaining the actual thickness but letting 
the density approach infinity. in comparison with the value of D 2 , D 2 _ Q 
is negligible. However, for the sake of comparison with exposure rate 
for other configurations having walls of mass thickness approaching 
infinity, one must consider ir. the limit the variation, with respect o 
direction, of the differential contribution of radiation entering the 
detector. The situation is one which is extremely difficult to establish 
experimentally, but there are theoretical means of obtaining radiation 
data at very great penetration depths; and, since the angular and energy 
distribution of the radiation field fends to seek an approximate equilib¬ 
rium at great depths, there is reason to expect that a description of 
can be obtained in such a way. The primary information needed is 

2-eo 
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the directional distribution of radiation emitted from the wall, or - 
looking at the matter from the detector point of view - the directional 
distribution of differential exposure contributions at the detector. 

NBS-76 implicitly assumes that the directional density function depicting 
the radiation emitted from the walls is separable into components depen¬ 
dent on vertical (polar) and horizontal (azimuthal) angles, respectively. 
Thus, this function is assumed to have the form: 

F(cos8,cp) = p(cp) q(cos0) , (6) 

where 0 and cp are the polar and azimuthal angles, respectively rela¬ 
tive to a frame of reference with a vertical polar axis. If the azimuthal 
angle is measured relative to a line normal to a wall, then NBS-76 gives 
as an expression for p(cp) , based upon an assumed cosine distribution, 



It is to be noted that the adequacy of the assumption as to separability 
of the function F(9,cp) has never been studied; and, even if it is 
approximately separable, the accuracy of the cosine assumption for p(cp) 
has not been verified. 

Situation 3-0 involves a circular structure of zero mass thickness; 
and thus D 3 _ Q is the exposure rate above an infinite contaminated plane, 
on the axis of a decontaminated circle which subtends the same solid angle 
at the detector position as does the rectangular base of the structure 
under consideration. A general formula for the exposure rate can be 
expressed as: 
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D 


3-0 



■t(H,cos0) dfi 


( 10 ) 


where the integration is over the solid angle subtended at the detector 
by the walls. Through the cylindrical symmetry of the situation, this 
becomes: 


0 

D 3-0 = j 


COS0 

f 1 

<6(H,cos0) d(cos0) + I {.(H,cos0) d(cos0) (11) 


-cos© 





-cos© 




-t(H,cos@) d(cos0) 


( 12 ) 


This formulation is not exact in principle, since the {.-function is 
for a completely contaminated plane and a collimated detector, while 
this case is for a partially cleared plane and a completely uncollimated 
detector. The degree of approximation is quite close, however, in most 
practical cases. 


Situation 3 -°q has a rectangular plan similar to that of the structure 
itself, but it is considered of infinite extent in the vertical direction. 
The mass thickness approaches infinity. The statements relative to Situa¬ 
tion 2-oo apply in general here also. Under the same assumptions, the 
exposure rate is found to be: 


j. 

L 


D = 2 J q(cos0) d(cos0) 
3-oo 


/ 


TT 

2 


= 2(sin9 coscp 

8 


-i 


coscp dcp + 

0 *"0 
q(cos0) d(cos0) 


<P« 


cos(p dcp 


(13) 


(14) 


Situation 3 ! -oe is similar to Situation 3-», except that the plan is 


circular instead of rectangular. The circle has the same area as that 
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specified for Situation 3-0. NBS-76 makes the same assumption of sepa¬ 
rability of vertical and horizontal dependence, and in this case an 
additional assumption is made that the vertical dependence is similar 
to that of the skyshine portion of the ^-function, that is: 

q(cos0) = k • £(3*,-]cos01) , (15) 

where k is an appropriate constant. The dependence on the negative 
of the absolute value of cos0 shews that the contribution from above 
and below the detector plane is believed to be symmetric at deep pene¬ 
trations. (One should note that here, just as in Situation 3-oo, the 
wall is effectively of infinite depth as well as infinite height.) For 
this case any line from the detector is perpendicular to the wall, and 
thus from Equation (7), p(cp) equals 1/2 at the detector position. 
Then 


2n 

D s | 4. 

S' -00 I 2 

J o 

0 

r 


2k | -t(3* , cosS) d(cos9) 


-1 


= 2nk 1 -M3', cos0) d(cos0) 

-1 

= 2nk S(3') 


(16) 


(17) 

(18) 


where S is a function defined in NBS-42 (Equation 27.3^. 


Situation 4-0 differs from Situation 3-0 only in that the cleared 
area within the "phantom'' structure is rectangular rather than circular, 
oeing similar in plan to the structure under consideration. A general 
formula for the exposure rate is given by 


D 


4-0 



£(H,cos0) dQ 


(19) 


where the integration is over the solid angle subtended at the detector 

(15 ) 

by the walls of the structure. The technology as established assumes 
the identity of and D 4 _ q , although NBS-76 notes that in principle 

this is not true. Calculation of D 4 _ Q can be made, although it is 
tedious because the limits of integration are not simple functions of the 
polar coordinates which are suitable for the ^-function. 
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Situation 4 -qq is the situation for a finite-sized building of rectan¬ 
gular plan. It is the same shape as the structure under investigation, 
but it has walls of infinite density. The angular distribution of detec¬ 
tor response in such case has not been previously defined, but in principle 
it can be obtained by theoretical means, if one can show that deep pene¬ 
trating radiation approaches an equilibrium directional distribution. In 
terms of the function F(cos,cp) introduced in the discussion of Situa¬ 
tion 2-oo t the exposure rate is simply given by: 


D 


4-00 



F(cos9,cp) dQ 


( 20 ) 


where the integration is over all solid angle subtended by the walls at 
the detector. 


Situation 4 '-qq is similar to that in Situation 3'-oo , except that 
the upper and lower limits on the wall are finite. Thus we would have, 
according to the distribution assumed by NBS-76: 


D = rrk 

4 "*oo 


{/ 

-C" 

0 

I l(! 


t(3 ', cos0) d(cos0) 
COS0. 


3',cos0) d(cos0) 


-cos© 


M- 


( 21 ) 


= TTk 


G(l-cos0.) + G(l-cos0 ) 
l U 


( 22 ) 


where G is defined in NBS-76 (Section 4.4, line 25). The validity of 
this equation is based on several assumptions as noted for Situations 
3-oo and 3'-» . 

Situation 4 is the real structure, of finite dimensions, actual wall 
density, and rectangular shape. The exposure rate is expressible as an 
integral over all directions subtended by the walls at the detector of 
the directionally differential detector response function, but this 
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definition, though obvious, is not directly useful. A direct and exact 
calculation is not easily done; and it would require an elaborate Monte 
Carlo type of computation not suited for engineering purposes. 


B.3 Discussion and Conclusions 

Since R is defined as D^/D^ , the value of R may be established 
in terms of the exposure rates for the various situations already dis¬ 
cussed, according to the following analysis: 



One may consider that the ratio D /D , which involves true wall mass 

thicknesses, can be interpolated between the corresponding ratios for 

infinite and zero mass wall thicknesses, respectively, that is , 

D, /D n and D /D. . . The interpolation factor is called S , 
4-00 2-oo 4-0 2-0 w 

and thus 


D 


4 


D 


2 




NBS-76 assumes that 


S = 
w 


B (X)-l 
P 

B (X) 

P 


(24) 


(25) 


where B is the buildup factor for cobalt-60 normally incident on con- 
P 

Crete (or similar material). Such assumption involves two separate hypo¬ 
theses: (1) the only significant parameter affecting the value of 

is the actual mass thickness of the structural wails; (2) the dependence 
is given by the ratio of scattered to uncollided cobalt-60 radiation which 
can penetrate the walls. 


Thus, one obtains, following this approach: 
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If the assumption is made that 



D 

D 


4-00 

3~00 


4'-00 



(28) 


which involves some degree of approximation, however close, then Equation 
(1) follows, through the establishment of the following definitions: 


B 

w 


G 

s 



E 

s 


E 


d 



4'—CO 



00 


D 


3-0 


D 

D 

D 

D 

D 


2-0 

3~oo 

2 - oo 

4-0 

3- 0 


(29) 

(30) 

(31) 

(32) 

(33) 


The mathematical descriptions of these concepts follow from the 

definitions. Since D, is unity, B is obtained by the same formulation 

1 w 

as D (see Equation 3 and 4 and associated text). G becomes, from 

A S 

Equations (18) and (22), 


G 

s 


G(i-cos0.) + G(l-cos0 ) 

L ____ u —I 

2 S(3') 


(34) 




5 



l-cos0 ) 




G(l-cos0^) 




(35) 


since S(3') is approximately equal to 0.1 . From Equations (5) and 
(12), one obtains: 
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G 


d 



£(H,cos0) 


L(H) 


d(cos0) 


(36) 


From Equations (9) and (14), one derives: 

E = sincp + cosep . (37) 

5 3 3 


From Equations (10) and (19), we see: 


E 


d 



£(H,cos0) dn 



<L(H,cos0) df) 


(38) 


where the integrations are over the solid angle subtended at the detector 

by the walls for the rectangular and equivalent cylindrical situations, 

as indicated. NBS-76 (Section 4.5) suggests that this ratio can be obtained 

(16 ) 

from the work of Hubbell et al . . Although the situations considered 

by Hubbell differed slightly from those involved here, the results so 

obtained would probably be a close approximation. Present engineering 
(15) 

practice , however, is to Tssume this value to be unity. 

It is to be noted that the concepts involving essentially infinitely 

thick structures, G and E , involve situations which are probably 

s s 

impossible to reproduce experimentally. Independent determinations of 
these so-called "non-wall scattered" factors^ 5 \ or even their product, 
by experimental means seem out of the question. Indirect experimental 
determination of these factors by measuring the overall reduction factor 
cannot be achieved unless all the other factors in Equation (1) are 
accurately known and the various assumptions inherent in the development 
of that formula have been justified. 


The important assumptions and approximations which are subject to 
verification or possible improvement are: 

(a) That deeply penetrating radiation can be split into factors 
depending separately upon vertical and azimuthal angles. 
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(b) That the azimuthal variation for deeply penetrating radiation 
follows a cosine distribution. 

(c) Tha; the vertical variation for deeply penetrating radiation 
is approximated by the skyshine portion of the /-function. 

(d) That the shape factor E J is always very near unity. 

d 

(e) That the interpolation factor S is dependent only on wall 

w 

mass thickness and that it is approximately given by the 
formula presented as Equation (25). 

(f) That G as determined on the basis of cylindrical structure 

s 

situations are valid for rectangular structures. 

One might ask at this point whether it is worth seriously considering 
abandoning the whole approach to ground contribution to structure pene¬ 
tration typified by Equation (1). However, it must be observed that the 
present formulation has been shown by many experiments to be rather suc¬ 
cessful thus far and is now deeply ingrained into the minds of many 

fallout shelter analysts. The edifice thus appears to require bolstering 

* 

rather than replacement. 

B.4 Recommendations 

It is recommended that il further work on ground penetration into 
blockhouses be undertaken, it be directed along the following lines in 
the immediate future: 

(a) Make further efforts to pin down for all practical thick¬ 

nesses, either experimentally or theoretically or both, so as to give 
confidence in definitive results drawn from such work within an uncer¬ 
tainty on the order of perhaps 10?- . 

(b) Bo further work, either experimentally or theoretically or both, 

to obtain or improve values of G and E for phantom" buildings 

a (1 

(i. e., those of zero wall mass thickness), to the same order of accuracy 
as indicated above. 

V. Eisenhauer disagrees with this (personal communication). He would 
prefer a completely fresh approach, avoiding the interpolation procedure. 
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(c) Determine E and G (or their product if they are not sepa- 

s s 

rablv.’) by approved theoretical techniques, to within the same order of 
accuracy as above. 

(d) After the above factors are known with reasonable accuracy, 
resume experimental blockhouse work to look for the best values of S 

w 

and to determine the effect of various parameter changes on it. This 

could possibly be supplemented with Monte Carlo calculations us'ng the 

actual blockhouse configurations. If the other functions are known within 

an uncertainty of about 10% , one should be able to determine J to 

w 

within about 20% . 

It is not implied that previous experimental work in these matters 
has been wasted. However, it is believed that after items (a), (b), 
and (c) have been accomplished, previous work could be re-interpreted in 
the light of this approach, along with further work mentioned in sub- 
paragraph (d). 

Some justification for continuation of work on this basic, elemen¬ 
tary problem appears to be necessary, especially in view of the apparant 
adequacy of the approach for simple geometries, as indicated above. Such 
justification arises most strongly from the fact that a great many of 
the problems involving more complicated type of structures and source 
distribution* are solved by methods which spring from this basic method 
and assume its adequacy. It is well known that there are serious dis¬ 
crepancies between experimental values and engineering predictions in 
connection with many of these complex situations, such as the blockhouse 
basement problem t in-and-down"', the protection prediction for upper 
floors of multistory structures (’ in-and-up* ). the mutual shielding 
problem (’’finite source plane"'*, and others. In order to solve these 
difficult problems in fallout shielding methodology, tho logical frame¬ 
work of the rjore basic part of the technology must be sufficiently firm 
to assure that the troubles in more complex situations do not stem, even 
in part, from these shaky foundations. Therefore, more careful attention 
in research should be devo’ed to strengthening the fundamentals of the 
present fallout shielding technology. 
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